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Abstract 

The dynamics of phase transitions plays a crucial role in the so-called interface 
between high energy particle physics and cosmology. Many of the interesting re¬ 
sults generated during the last fifteen years or so rely on simplified assumptions 
concerning the complex mechanisms typical of nonequilibrium field theories. In 
particular, whenever first order phase transitions are invoked, the metastable 
background is assumed to be sufficiently smooth to justify the use of homoge¬ 
neous nucleation theory in the computation of nucleation rates of critical bubbles. 
In this talk I present tlie results of numerical simulations which were designed to 
quantify “smoothness”; that is, how the contribution from nonperturbative sub- 
critical fluctuations may spoil the homogeneity assumption of nucleation theory. 
I then show how the numerical results can be understood quantitatively in terms 
of a simple analytical model of subcritical thermal fluctuations. Encouraged by 
the success of the model in matching the numerical results, I apply it to the 
standard model electroweak phase transition. 
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I. INTRODUCTION: HOMOGENEOUS NUCLEATION AND ITS 


ASSUMPTIONS 


The fact that the gauge symmetries describing particle interactions can be restored at 
high enough temperatures has led, during the past 15 years or so, to an active research 
program on the possible implications that this symmetry restoration might have had to the 
physics of the very early Universe. One of the most interesting and popular possibilities 
is that during its expansion the Universe underwent a series of phase transitions, as some 
higher symmetry group was successively broken into products of smaller groups, up to the 
present standard model described by the product SU{3)c <8) SU{2)l ® U{1)y- Most models 
of inflation and the formation of topological (and nontopological) defects are well-known 
consequences of taking the existence of cosmological phase transitions seriously 

One, but certainly not the only, motivation of the works addressed in this talk comes 
from the possibility that the baryon asymmetry of the Universe could have been dynamically 
generated during a first order electroweak phase transition p| . As is by now clear, a realistic 
calculation of the net baryon number produced during the transition is a formidable chal¬ 
lenge. We probably must invoke physics beyond the standard model (an exciting prospect 
for most people) 0, push perturbation theory to its limits (and beyond, due to the nonper- 
turbative nature of magnetic plasma masses that regulate the perturbative expansion in the 
symmetric phase), and we must deal with nonequilibrium aspects of the phase transition. 
Here I will focus on the latter problem, as it seems to me to be the least discussed of the 
pillars on which most baryon number calculations are built upon. To be more specihc, it 
is possible to separate the nonequilibrium aspects of the phase transition into two main 
subdivisions. If the transition proceeds by bubble nucleation, we can study the propagation 
of bubbles in the hot plasma and the transport properties of particles through the bubble 
wall. A considerable amount of work has been devoted to this issue, and the reader can 
consult the works of Ref. ^ for details. These works assume that homogeneous nucleation 
theory is adequate to investigate the evolution of the phase transition, at least for the range 
of parameters of interest in the particular model being used to generate the baryon asymme¬ 
try. This brings us to the second important aspect of the nonequilibrium dynamics of first 
order phase transitions, namely the validity of homogeneous nucleation theory to describe 
the approach to eqnilibrium. This is the issue addressed in this talk. 

Nucleation theory is a well-studied, but far from exhausted, subject. Since the pioneering 
work of Becker and Doring on the nucleation of droplets in supercooled vapor , the study 
of first order phase transitions has been of interest to investigators in several fields, from 
meteorology and materials science to quantum field theory and cosmology. Phenomenologi¬ 
cal field theories were developed by Cahn and Hilliard and by hanger in the context of 
coarse-grained time-dependent Ginzburg-Landau models, in which an expression for the de¬ 
cay rate per unit volume was obtained by assuming a steady-state probability current flowing 
through the saddle-point of the free-energy functional The application of metastable 

decay to quantum held theory was initiated by Voloshin, Kobzarev, and Okun [Q , and soon 

The generalization of 
and has been 


after put onto hrmer theoretical ground by Coleman and Callan 
these results for hnite temperature held theory was hrst studied by Linde |]n 
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the focus of much recent attention 12 


The crucial ingredient in the evaluation of the decay rate is the computation of the 
imaginary part of the free energy. As shown by hanger [|^ , the decay rate TZ is proportional 
to the imaginary part of the free energy JF, 


= 


E. 


ttT 


-ImjF , 


( 1 ) 


where is the negative eigenvalue related to metastability, which depends on nonequi- 
librium aspects of the dynamics, such as the growth rate of the critical bubble. Since 
T = —TlnZ, where Z is the partition function, the computation for the rate boils down to 
the evaluation of the partition function for the system comprised of critical bubbles of the 
lower energy phase inside the metastable phase. 

If we imagine the space of all possible held conhgurations for a given model, there will be 
different paths to go from the metastable to the ground state. We can think of the two states 
as being separated by a hill of a given “height”. The energy barrier for the decay is then 
related to the height of this hill. At the top of the hill, only one direction leads down to the 
ground state, the unstable direction. Fluctuations about this direction will grow, with rate 
given by the negative eigenvalue which appears in the above formula. All other directions 
are positively curved, and fluctuations about them give rise to positive eigenvalues which 
do not contribute to the decay rate. The path which will cost less energy is the one which 
will dominate the partition function, the so-called critical bubble or bounce. It is simply the 
held conhguration that interpolates between the two stable points in the energy landscape, 
the metastable and ground state. The energy barrier for the decay is the energy of this 
particular held conhguration. 

For a dilute gas of bubbles only, the partition function for several bubbles is given by 
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Z ~ Z{^f) + Z{ipf) 
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where </?/ is the metastable vacuum held conhguration and (pb is the bubble conhguration, 
the bounce solution to the 0(3)-symmetric Euclidean equation of motion. We must evaluate 
the partition functions above. This is done by the saddle-point method, expanding the scalar 
held (/)(x, r), such that (/)(x, r) —*• -|- C(x, r) for Z{ipf), and 0(x, r) —>• (pb{^) + ? 7 (x, r) for 

Z{Lpb)i where C(x, r) and ? 7 (x, r) are small huctuations about equilibrium. 

It is crucial to note that the saddle-point, or Gaussian, method only gives good results if 
indeed the huctuations about equilibrium are sufficiently small that nonlinear terms in the 
helds can be neglected. Even though the method sums over all amplitude huctuations, it does 
so by assuming that the functional integral is well approximated by truncating the expansion 
of the action to second order. The efficiency of the method relies on the fact that higher 
amplitudes will be suppressed fast enough that their contribution to the partition function 
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will be negligible. One can visualize this by comparing a sharp parabolic curve with a flatter 
one with minimum at xq, and investigating when / will be well approximated by 

writing /(x) ~ /(xq) + |(x — xqY f"{xo). For a sharp curve, larger amplitude fluctuations 
will be strongly suppressed and thus give a negligible contribution to the integral over all 
amplitudes. Clearly, this will not be the case for flatter curves. 

Skipping details [|^|, using the saddle-point method one obtains for the ratio of partition 
functions, |[g}. 


^ (^b) saddle—point 

Z{^f) 


det(-aj; V”{iph))f3 
det(-nE V"{(pf))p 


.-A5 


( 3 ) 


where [det(M)^] 2 = / F)r 7 exp|— /(f dr f d^x\r][M]rj^ and /S.S = SE^^h) — is the 

difference between the Euclidean actions for the held configurations and ipf. [Note that 
Se{^), and hence AS, does not include any temperature corrections. It would if one had 
summed over other fields coupled to ip.] Thus, the free energy of the system is. 


T = -T 


det(-ai; -F V''{pi,))p 
det(-ni? -h V''{pf))fs 


,-AS 


( 4 ) 


Let me stress again the assumptions that go into computing the free energy. First, that 
the partition function is given by Eq. ^ within the dilute gas approximation, and second, 
that the partition function is evaluated approximately by assuming small fluctuations about 
the homogeneous metastable state pf. It is clear that for situations in which there are 
large amplitude fluctuations about the metastable equilibrium state the above formula must 
break down. Thus the breakdown of the expression for the rate is intimately connected 
with the question of how well-localized the system is about the metastable state as the 
temperature drops below the critical temperature T^. Homogeneous nucleation, as its name 
already states, is only accurate when the metastable state is sufficiently homogeneous. In 
the presence of inhomogeneities, there is no reason to expect that the decay rate formula 
will apply. The question then is to quantify when does it breakdown. 

In order to study the properties of the metastable background in the presence of ther¬ 
mal fluctuations, I have recently investigated the nonequilibrium dynamics of a (2-1-1)- 
dimensional model of a scalar held coupled to a thermal bath |^. The nonlinear interactions 
were chosen to reflect the gross properties of the electroweak effective potential, although 
the model only deals with a real scalar held. These results were then extended to a (3-1-1)- 
dimensional simulation, in work done in collaboration with Julian Borrill |T^. The results 


for two and three dimensions were qualitatively similar, showing how the amount of phase 
mixing is related to the parameters of the model. In particular, for the three-dimensional 
study we varied the strength of the scalar field’s quartic coupling, which for the electroweak 
model is related to the Higgs boson mass. In the next Section I review the results of this 
work. In Section 3 I present an improved treatment of the kinetics of subcritical bubbles 
which I then use to understand the results of the numerical simulations. The agreement 
between numerical experiment and analytical model is extremely satisfactory. The simplic¬ 
ity of the kinetic approach based on the dynamics of subcritical thermal fluctuations makes 
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it a useful tool in studying the properties of the thermal background. I then briefly apply 
the subcritical bubbles method to the simple 1-loop electroweak potential, showing that for 
Higgs masses above 60 GeV, the amount of phase mixing is such that homogeneous nucle- 
ation should not be used to study the dynamics of the transition. I conclude in Section 4 
with a brief overview of the results and possible future work. 


II. NUMERICAL SIMULATIONS OF THERMAL PHASE MIXING 


The homogeneous part of the free energy density is written as 


r) = ^ (r= - r|) 4 ,^ - + G 


(5) 


This choice intentionally resembles the electroweak effective potential to some order in per¬ 
turbation theory, although here 0(x, t) is a real scalar held, as opposed to the magnitude of 
the Higgs held. The goal is to explore the possible dynamics of a model described by the 
above free-energy density, generalizing the results obtained in Ref. to (3-|-l)-dimensions. 
The analogy with the electroweak model is suggestive but not quantitative. 

Introducing dimensionless variables x = a^T 2 X, i = a^T 2 t, X = a~iT 2 ^ 4 >, and 9 = T/T 2 , 
the Hamiltonian is. 


H[X] 


1 

0 


(fx 


i I VA- P +i (9= - 1 ) A= - 


“ex’ + -X* 

3 4 


( 6 ) 


3 ~ 1 

where a = a~^a, and A = a~ 2 \ (henceforth we drop the tildes). For temperatures above 
61 = [1 — Q!^/ 4 A )“2 there is only one minimum at X = 0. At 6 = 61 an inhection point 
appears at Xinf = a6i/2X. Below 61 the inhection point separates into a maximum and a 
minimum given by X± = ^ 1 ± y 1 — 4A (1 — 1/0^) . At the critical temperature 9c = 


(1 — 2a^/9A)“5 the two minima, at Xq = 0 and X_|_ are degenerate. Below 9c the minimum 
at X_|_ becomes the global minimum and the Xg-phase becomes metastable. Finally, at 0 = 1 
the barrier between the two phases at X_ disappears. 

The coupling with the thermal bath will be modelled by a Markovian Langevin equation 
which, in terms of the dimensionless variables dehned above, is 


d^X 

dt'^ 


dX dU(X,9) , 


(7) 


where rj is the dimensionless viscosity coefficient, and ^ the dimensionless stochastic noise 
with vanishing mean, related to r] by the huctuation-dissipation theorem. 


(^(x, t)^(x', t')) = 2T]95{t - - x') 


( 8 ) 


A few comments are in order concerning our choice of equation. It is clear that we are 
assuming that X (x, t) represents the long-wavelength modes of the scalar held. Whenever 
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one discretizes a continuum system there is an implicit coarse-graining scale built in. We 
encapsulate information about the shorter-wavelength modes, which have faster relaxation 
time-scales, in the dissipation and noise terms. In principle it should be possible to derive an 
effective Langevin-like equation for the slow modes by integrating out the fast modes from 
the effective action. This is a complicated problem, and progress has been slow. Recent work 
indicates that one should expect departures from the Langevin equation written above [ITB 


although details are sensitive to the particular model one starts with. For example, the noise 
may be colored (with more complicated correlation functions) and the coupling to the bath 
may be multiplicative, as opposed to the additive coupling chosen above. Here, we will adopt 
the above equation as a hrst step. We do not expect that the nature of the noise will change 
the hnal equilibrium properties of the system, but mostly the relevant relaxation time-scales. 
Furthermore, Gleiser and Ramos showed that for high enough temperatures the noise does 
become white |^. Since the physical results here are related to the hnal equilibrium state 


of the system, we believe that they will not be affected by more complicated representations 
of the coupling of the held to the thermal bath. However, a more thorough examination of 
this question deserves further study. 

A related topic is the choice of coarse-graining scale, which is embedded in the lattice 
spacing used in the simulations. It is well-known that any classical held theory in more 
than one spatial dimension is ultra-violet divergent, and that the lattice spacing serves as an 
ultra-violet cutoh. This being the case, one should be careful when mapping from the lattice 
to the continuum theory. If one is to probe physics at shorter wavelengths, renormalization 
counterterms should be included in the lattice formulation so that a proper continuum limit 
is obtained on the lattice within the validity of perturbation theory. This point has been 
emphasized in Ref. []^, where a (2-|-l)-dimensional study of nucleation was performed for 
a temperature-independent potential. Renormalization counterterms (of order 6 In 5a; for 
lattice spacing 6x) for a particular renormalization prescription were obtained, and the 
results shown to be lattice-space independent. 

Here, due to the temperature dependence of the potential, the renormalization prescrip¬ 
tion of Ref. [Q does not work. Instead, we will use 5a; = 1 throughout this work. It turns 
out that for all cases studied the mean-held correlation length = l/"(Xo, 9c) is suffi¬ 
ciently larger than unity to justify this choice. Modes with shorter wavelengths are coupled 
through the noise into the dynamics of the longer wavelength modes, as described by Eq. 
above. But it should be stressed that different choices of lattice spacing imply in different 
effective coarse-graining scales, and hence different results. In other words, although phase 
mixing will always be present to some degree, the exact quantitative results will depend on 
the coarse-graining scale. 

I will skip details of how we implemented this simulation in a parallel machine, as well 
as the tests to make sure our results were independent of simulation parameters such as 
lattice length, time and spatial steps, and random number generator. The reader is referred 
to Ref. |T^ for information on these issues. I will now proceed by describing the numerical 


experiments and our results. 

As pointed out in the Introduction, the question of whether homogeneous nucleation the¬ 
ory is trustworthy to describe a hrst order phase transition boils down to how well localized 
in the metastable state the system is, as the temperature drops below the critical tem- 
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perature. In the jargon of condensed matter physics, nucleation should work for quenches 
deep into the metastable branch of the coexistence curve. In order to address this question, 

HI 


following the procedure of Ref. 


we will study the behavior of the system at the crit¬ 


ical temperature, when the two minima are degenerate. The reason for this choice follows 
naturally from the fact that we are interested on the way by which the system approaches 
equilibrium as the temperature drops below Tc. The detailed dynamics will depend on the 
relative fraction of the total volume occupied by each phase; if at T^. the system is well 
localized about the X = 0 minimum, as the temperature drops the transition may evolve by 
nucleation and subsequent percolation of bubbles larger than a critical size. If, on the other 
hand, considerable phase-mixing occurs already at Tc, we expect the transition to evolve 
by domain coarsening, with the domains of the phase eventually permeating the whole 
volume. 

Let us call the two phases the 0-phase and the -I— phase, corresponding to the local 
equilibrium values X = Xq = 0, and X = X^, respectively. We can quantify the phase 
distribution of the system as it evolves according to Eq. |^, by measuring the fraction of the 
total volume in each phase. This is done by simply counting the total volume of the system 
at the left of the potential barrier’s maximum height (he., X < = X ma ^), corresponding 

to the 0-phase. Dividing by the total volume, we obtain the fraction of the system in the 
0-phase, fo(t), such that 


fo(t) + f+(t) — 1 , (9) 

where, of course, f+{t) corresponds to the fractional volume in the -I— phase. A further 
measure of any conhguration is given by the volume-averaged order parameter, (X)(t) = 
V~^ f dV X(t). A localized configuration > 0.5) then corresponds to (X)eq < X^ax, and 
a fully phase-mixed conhguration (/q*^ ~ 0.5) to (X)eq = Xmax, where the super(sub)-script 
‘eq’ refers to hnal ensemble-averaged equilibrium values of fo(t) and (X)(t), respectively. 
[Recall that in the presence of noise we must take an ensemble average of the physical 
quantities of interest in order to get smooth, sensible results. In practice, this amounts to 
changing the seed of the random number generator for each run, and averaging the results 
over many runs.] 

We prepare the system so that initially it is well localized in the 0-phase, with /o(0) = 1 
and (X)(0) = 0. These initial conditions are clearly the most natural choice for the problem 
at hand. If one has cosmology in mind, it is quite possible that as the system slowly cools 
down (we are not interested in phase transitions close to the Planck scale), fluctuations 
from the high temperature phase X = 0 to the X_|_ phase are already occurring before Tc is 
reached. (In this case, our arguments are even stronger.) However, we will adopt the best- 
case scenario for homogeneous nucleation to work, in which the system managed to reach 
the X = 0 phase homogeneously, so that the initial state is a thermal state with mean at Xg. 
If one has more concrete applications in mind, we can assume that we quenched the system 
to its critical temperature, making sure that the order parameter remains localized about 
the high-temperature phase. Since thermalization happens very fast in the simulations, the 
exact point by point initial conditions should not be important, and we can view the first few 
time steps as generating an initial thermal distribution with /o(0) ~ 1 and (X) ~ 0, so that 
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the average kinetic energy per lattice point satisfies the equipartition theorem, = |T. 
For simplicity we take X = 0,X = 0 everywhere initially. 

There are two parameters controlling the strength of the transition, a and A. In the 
previous (2+l)-dimensional work, a was chosen to vary while A was kept fixed. It is really 
immaterial which parameter is held fixed, or if both are made to vary, but in order to keep 
closer to the spirit of the electroweak model we will fix a and let A vary. As is well-known, A 
is related to the Higgs mass, while a is related to the gauge-boson masses 0 . The connection 
with the electroweak model is straightforward. If we consider as an example the unimproved 
one-loop approximation, the effective potential is |^, 

Hew(0, T) = D (t^ - T|) 02 _ eT<P^ + 1at0", (10) 

where D and E are given by iA = [^{Mw/aY + 3{Mz/aY + 6(Me/ct) 2] /24 ~ 0.17 and E = 
[6{Mw/(tY + 3{Mz/aY] /12n ~ 0.097, for Mw = 80.6 GeV, Mz = 91.2 GeV, Mt = 174 
GeV 0 , and a = 246 GeV. T 2 is given by. 


T 2 = ^(M2 - 8Ha2)/4ZA , (11) 

where the physical Higgs mass is given in terms of the 1-loop corrected A as Mjj = 
(2A -I- 12i?) (j2, with B = (6M^ -|- 3M| — 12M|.) /bdvr^cr'^, and the temperature-corrected 


Higgs self-coupling is. 


Xt — A — 


Ibvr^ 




B 




B 
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In 


Mp 


a 


In (mI/cfT^) 


( 12 ) 


where the sum is performed over bosons and fermions (in our case only the top quark) with 
their respective degrees of freedom gB{F), and Inc^ = 5.41 and Incp = 2.64. 

Thus, the correspondence with our (dimensionless) parameters is 


3E 

(2D)T 


0.065 , and A 


Xf 

TFT 


1.72At 


(13) 


Once this is established, the numerical experiment proceeds as follows: i) Ghoose a = 0.065; 
ii) Prepare the system in the initial state described above, and measure the value of /o(f) 
and {X)(t) for several values of A, as the system evolves according to Eq. |^. I can now 
present the results of the simulations. 

Based on the above discussion, we choose lattice length L = 48, lattice spacing 6x = 1, 
time step 6t = 0.1, and a = 0.065 in all simulations. [For details as to why this choice of 
parameters is sensible, see Ref. [|^.] The experiment then consists in measuring the fraction 


of the volume in the 0-phase as a function of time for several values of A. 

In Fig. 1 we show the evolution of the ensemble-averaged fraction fo{t) for several 
values of A. It is clear that for small enough values of A the system remains well-localized 
in the 0-phase with fY ~ 1, while for larger values the two phases become completely 
mixed, with /q^ —> 0.5. Remarkably, the transition region between the two regimes is quite 













narrow, centered around A ~ 0.025. This can be seen from Fig. 2 where we show fo{t) for 
A = 0.024, 0.025, and 0.026. [The curves are noisier due to the fact that we must run for 
longer times in order to approach the equilibrium values of fo{t), being thus constrained to 
perform an ensemble average with fewer runs.] 


Figure 1: The approach to equilibrium for several values of A. 


Figure 2: Fitting fo{0c) to a power law at large times for A = 0.025. 


Note that for A = 0.026, /g’’ ~ 0.5, while for A = 0.024, ~ 0.72. There is a 

pronounced change in the behavior of the system for A ~ 0.025. Furthermore, we hnd that 
the numerical curves can be htted at all times by a stretched exponential, 

foit) = (1 - fP) exp [-(t/req)'"] + /o^ , (14) 

where /g*^ is the dual equilibrium fraction and Xeq is the equilibration time-scale. In Table 1 
we list a and Teq for several values of A. Note that for A = 0.025 the fit is obtained at late 
times by a power law, (smooth curve in Fig. 5) 

foit) |a=0.025 oc t ^ , (15) 


with k = 0.10(±0.02). This slowing down of the approach to equilibrium is typical of 
systems in the neighborhood of a second order phase transition, being known as ‘critical 
slowing down’ [^. This behavior is suggestive of a “phase transition” between two possible 
regimes for the system, one in which the system is well-localized in the 0-phase, and the 
other in which there is a complete mixing between the two phases. Let us call these two 
regimes the ‘smooth’ and the ‘mixed’ regimes, respectively. [In Refs, fl^ and [|T^] I used 
the terms “strong” and “weak”, respectively. The change is not to induce confusion, but 
simply to emphasize that the difference between the two regimes is in the amount of mixing 
present in the background. Strong and weak refer to the strength of the associated phase 
transition, although I consider the term weak hrst order transition a misnomer; the question 
is nucleation vs. spinodal decomposition and not weak vs. strong.] Before we explore this 
idea any further, I note that the final equilibrium fractions are insensitive to the viscosity 
parameter r], which reflects the coupling of the system to the thermal bath. For details. 
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I refer the reader to Ref. [15|. This is precisely what one expects, as the coupling to the 


bath should not influence the final equilibrium properties of the system, but only how fast 
it equilibrates. 


Table 1; The values of the equilibration time-scales and the exponents for the exponen¬ 
tial £t of Eq. 0 for several values of A. Uncertainties are in the last digit. 


Armed with these results, and invoking also the results in (2 -|- l)-dimensions 
dehne the equilibrium fractional population difference 
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we 


AFeq (0c) =/o" -. (16) 

In Fig. 3 we show the behavior of ATeq as a function of A. There is a clear qualitative anal¬ 
ogy between the behavior of ATeq as a function of A and the behavior of the magnetization 
as a function of temperature in Ising models. Here, the order parameter is the equilibrium 
fractional population difference and the control parameter is A. Ac — 0.025 is the critical 
value for the parameter A, which determines the degree of mixing of the system at T^.. 


Figure 3: The fractional equilibrium population difference ATeq as a function of A. 


We stress that the idea here is to probe the assumption of localization within the 0-phase 
as the system cools to Tc. Our results show that if the time-scales for cooling are slower 
than the equilibration time-scales of the system, for A > Ac there will be considerable phase 
mixing before the temperature drops below T^. This result can be made quite transparent 
by comparing the equilibrium value of the volume-averaged held (A)eq, and the location of 
both the inhection point and the maximum of the potential with varying A. As can be seen 
from Fig. 4, the narrow transition region is clearly delimited by 

Ainf < (W)eq < Amax , (17) 

where Xinf and are the inhection point and the maximum of the potential barrier, 

respectively. Note that for A > 0.026, = 0.5 and (X)eq = Xmax- 
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Figure 4: Comparison between volume-averaged field and location of the inflection point 
Xinf and top of the barrier X^ax, as a function of A. 


Recalling the information from Fig. 3, we conclude that there is a clear distinction 
between the ‘smooth’ and ‘mixed’ regimes. The fact that the mean-field potential used in 
the simulation has a barrier between the two phases does not imply that we can assume that 
the background is sufficiently smooth in order for homogeneous nucleation to be applicable. 
Thermal fluctuations, not included in the coarse-grained potential, change the character of 
the transition by promoting an efficient mixing between the two stable phases. In other 
words, the effective potential describing the true behavior of the system turns into a single¬ 
well potential with minimum at X^ax for A > Ac, even though the mean-field potential has a 
barrier between the two phases. Similar behavior exists in ferromagnets, where it is known 
that the true Curie temperature is below that predicted by mean-field theory. In order to 
improve on the mean-field results, one uses large N or ^-expansion methods, which sum 
over a larger class of diagrams incorporating leading infrared effects. It is possible that our 
results could be understood in part by using e-expansion methods. However, in the next 
Section I will propose a different approach to studying phase mixing, based on a kinetic 
equation for subcritical nonperturbative fluctuations, which, I believe, are responsible for 
the mixing. This is a new version of the subcritical bubbles method proposed by Gleiser, 
Kolb, and Watkins in Ref. [^. It will be shown that this method provides a simple and 
powerful way to study the background. 


III. MODELING PHASE MIXING WITH SUBCRITICAL BUBBLES 


As was stressed before, the computation of decay rates based on homogeneous nucleation 
theory assumes a smooth metastable background over which critical bubbles of the lower free 
energy phase will appear, grow and coalesce, as the phase transition evolves. However, as the 
results from the numerical simulation indicate, the assumption of smoothness is not always 
valid. To the skeptical reader, I point out that several condensed matter experiments indicate 
that homogeneous nucleation fails to describe the transition when the nucleation barrier 
(AS'/T) becomes too small. Furthermore, the agreement between theory and experiment 
has a long and problematic history [^. Homogeneous nucleation has to be used with care, 
in a case by case basis. It is thus important to obtain a criterion which will allow us to 
specify the conditions when homogeneous nucleation theory does apply. Below I describe a 
simple method which does exactly that. 

The basic idea is that in a hot system, not only small but also large amplitude fluctua¬ 
tions away from equilibrium will, in principle, be present. Small amplitude fluctuations are 
perturbatively incorporated in the evaluation of the finite temperature effective potential. 
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following well-known procedures. Large amplitude fluctuations probing the nonlinearities of 
the theory are not. Whenever they are important, the perturbative effective potential be¬ 
comes unreliable. In an ideal world, we should be able to sum over all amplitude fluctuations 
to obtain the exact partition function of the model, and thus compute the thermodynamic 
quantities of interest. However, we can only to this perturbatively, and will always miss 
information coming from the fluctuations not included in its evaluation. If large amplitude 
fluctuations are strongly suppressed, they will not contribute to the partition function, and 
we are in good shape. But when are they important? We can try to approach this question 
avoiding complicated issues related to the evaluation of path integrals beyond the Gaus¬ 
sian approximation by obtaining a kinetic equation which describes the fraction of volume 
populated by these large amplitude fluctuations. In order to keep the treatment simple, 
and thus easy to apply, several assumptions are made along the way, which I believe are 
quite sensible. In any case, the strength of the method is demonstrated when the results are 
compared with the numerical experiments described before. 

Large amplitude fluctuations away from equilibrium are modelled by Gaussian-prohle 
spherically-symmetric held conhgurations of a given size and amplitude. They can be 
thought of as being coreless bubbles. Keeping with the notation of the numerical exper¬ 
iment, huctuations away from the 0-phase, and into the 0-phase are written respectively 
as. 


Mr) = M , Mr) = </>c (l - e , (18) 

where R is the radial size of the conhguration, and (pc is the value of the amplitude at the 
bubble’s core, away from the 0-phase. In previous treatments (cf. Refs, 
was assumed that (pc = (p+i that is, that the conhguration interpolated between the two 
minima of the ehective potential, and that R = ^{T), where ^{T) = m{T)~^ is the mean- 
held correlation length. But in general, one should sum over all radii and amplitudes above 
given values which depend on the particular model under study. This will become clear as 
we go along. 

Dehne dn{R, (p, t) as the number density of bubbles of radius between R and R + dR at 
time t, with amplitudes (p > (pc between (p and (p -1- dcp. By choosing to sum over bubbles of 
amplitudes (pc and larger, we are ehectively describing the system as a “two-phase” system. 
For example, in the numerical simulation above it was assumed that the 0-phase was for am¬ 
plitudes (p < 0max, and the -I— phase was for amplitudes (p > 0max- Glearly, for a continuous 
system this division is artihcial. However, since the models we are interested in have two 
local minima of the free energy, this division becomes better justihed. Fluctuations with 
small enough amplitude about the minima are already summed over in the computation of 
the effective potential. It is the large amplitude ones which are of relevance here. To sim¬ 
plify the notation, from now on I will denote by “-1—phase” all fluctuations with amplitudes 
(p > (pc and larger. The choice of (pc is model-dependent, as will be clear when we apply this 
formalism to specihc examples. 

The fact that the bubbles shrink will be incorporated in the time dependence for the 


M and ra), it 
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radius jtS Here, I will only describe a somewhat simplified approach to the dynamics. More 
details are provided in the forthcoming paper by Gleiser, Heckler, and Kolb The results, 
however, are essentially identical. 

The net rate at which bubbles of a given radius and amplitude are created and destroyed 
is given by the kinetic equation. 


dn{R, t) 

Ft 



(19) 


Here, ro^+(i?) (r+^o(-R)) is fhe rate per unit volume for the thermal nucleation of a bubble 
of radius R of +-phase within the 0 -phase ( 0 -phase within the -I— phase). Ttn{R) — oT / 
is the (somewhat ad hoc) expression used for the thermal destruction rate, with a a constant 
related to the coupling to the thermal bath. Vo(+) is the volume available for nucleating 
bubbles of the -|-(0)-phase. Thus we can write, for the total volume of the system, Vt = 
Vo + h+, expressing the fact that the system has been “divided” into two available phases, 
related to the local minima of the free energy density. It is convenient to define the fraction 
of volume in the -I— phase, 7, as 


( 20 ) 

In order to compute 7 we must sum over all bubbles of different sizes, shapes, and 
amplitudes within the -I— phase, he., starting with > 0c- Clearly, we cannot compute 7 
exactly. But it turns out that a very good approximation is obtained by assuming that the 
bubbles are spherically symmetric, and with radii above a given minimum radius, i?min- The 
reason we claim that the approximation is good comes from comparing the results of this 
analytical approach with numerical simulations. The approximation starts to break down 
as the background becomes more and more mixed, and the morphology of the “bubbles” 
becomes increasingly more important, as well as other terms in the kinetic equation which 
were ignored. For example, there should be a term which accounts for bubble coalescence, 
which increases the value of 7 . This term becomes important when the density of bubbles 
is high enough for the probability of two or more of them coalescing to be non-negligible. 
As we will see, by this point the mixing is already so pronounced that we are justified in 
neglecting this additional complication to the kinetic equation. As a bonus, we will be able 
to solve it analytically. The expression for 7 is. 


^Of course, the amplitude 0 will also be time-dependent. However, its time-dependence is coupled 
to that of the radius, as recent studies have shown [24|. In order to describe the effect of shrinking 
on the population of bubbles it is sufficient to include only the time dependence of the radius. 
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r°° ( \ d^n 


0 min '^^ m in \ 3 / d(j)dR 


d(j)dR 


( 21 ) 


The attentive reader must have by now noticed that we have a coupled system of equa¬ 
tions; 7 , which appears in the rate equation for the number density n, depends on n itself. 
And, to make things even worse, they both depend on time. Approximations are in order, 
if we want to proceed any further along an analytical approach. The first thing to do is to 
look for the equilibrium solutions, obtained by setting dn/dt = 0 in the kinetic equation. In 
equilibrium, 7 will also be constant|. If wished, after finding the equilibrium solutions one 
can hnd the time-dependent solutions, as was done in Ref. . Here, we are only interested 
in the final equilibrium distribution of subcritical bubbles, as opposed to the approach to 
equilibrium. 

The first approximation is to take the shrinking velocity of the bubbles to be constant, 
dR/dt = —V. This is in general not the case (c/. Ref. [^), but it does encompass the fact 
that subcritical bubbles shrink into oblivion. The strength of the thermodynamic approach 
is that details of how the bubbles disappear are unimportant, only the time-scale playing 
a role. The second approximation is to assume that the rates for creation and destruction 
of subcritical fluctuations are Boltzmann suppressed, so that we can write them as T = 
^ 2 n 4 g-Fsc/T, ^ jg arbitrary constant of order unity, and Fsc{R,(l)c) is the cost in 

free energy to produce a conhguration of given radius R and core amplitude 0c- For the 
Gaussian ansatzwe are using, Fgc assumes the general form, F^c = aR + ^R^, where a = bcj)^ 
{b is a combination of tt’s and other numerical factors) and (3 depends on the particular 
potential used. In practice, the cubic term can usually be neglected, as the free energy of 
small (i? ~ subcritical bubbles is dominated by the gradient (linear) term. We chose to 
look at the system at the critical temperature Tc. For this temperature, the creation and 

destruction rates, ro^+ and T^_ ,0 are identical. Also, for Tc, the approximation of neglecting 

the cubic term is very good (in fact it is better and better the larger the bubble is) even for 
large bubbles, since for degenerate vacua there is no gain (or loss) of volume energy for large 
bubbles. Finally, we use that R+/Fr = 7 in the F+^q term, and that f+/Fr = n{R)V^{R) 
in the Ftd term. This latter expression assumes that the distribution of bubbles is sharply 
peaked around the smallest bubbles considered, so that the sum over all radii is fairly well 
approximated by one term. A more sophisticated approach is presented in Ref. We 

can then write the equilibrium rate equation as. 


nil 

^^dn(R)-cf(R) , 


( 22 ) 


where. 


^This doesn’t mean that thermal activity in or between the two phases is frozen; equilibrium is 
a statement of the average distribution of thermodynamical quantities. Locally, bubbles will be 
created and destroyed, but always in such a way that the average value of n and 7 are constant. 
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■Fsc/T 


d = aT/v, c = {1 — 2'j)AT'^/v, f{R) = e 


( 23 ) 


Integrating from -Rmin and imposing that n{R —>• cx)) = 0, the solution is easily found to 
be, 


n{R) 


^ -ai<t,,)R/T 

{d + a{(pc)/T) 


(24) 


Not surprisingly, the equilibrium number density of bubbles is Boltzmann suppressed. 
But we now must go back to 7 , which is buried in the dehnition of c. We can solve for 7 
perturbatively, by plugging the solution for n back into Eq. After a couple of fairly 
nasty integrals, we obtain. 


g (cr(0min)) -Rmin) ) 

-Rmin) 


(25) 


where. 


9 (®(0min)) Rmin) 


dvr /AT^\ 

^ V W 


(T\^ e r faRram\^ , „ai?min W , aRmm\ 

U) b + —) 


(26) 


A. Comparison with the numerical simulations 

We can now apply this formalism to any model we wish. The hrst obvious application 
is to compare 7 obtained from the numerical experiment with the value obtained from the 
kinetic approach. From the dehnition of the equilibrium fractional population difference, 
Eq. 0, 


AFeq (0e) = 1 - 27 . (27) 

Thus, it is straightforward to extract the value of 7 from the numerical simulations, as a 
function of A. Also, as we neglected the volume contribution to the free energy of subcritical 
bubbles, we have. 


3\/2 

F,, = a(0,)i?^i„ = ^7r3/2x!(0,)i?^i, , (28) 

O 

where, as you recall, X_ is the position of the maximum of the mean-held potential used 
in the simulations. So, we must sum over all amplitudes with X > X_, and all radii with 
7? > 1 (in dimensionless units), as we took the lattice spacing to be £ = 1. That is, we 
sum over all possible sizes, down to the minimum cut-oh size of the lattice used in the 
simulations. In practice, we simply substitute 0c = X_ and i?min = 1 in the expression 
for 7 . In Fig. 5, we compare the numerical results for 7 (dots) with the results from the 
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analytical integration of the kinetic eqnation. We took Ajv = 1, and the two cnrves are for 
ajv = 0 (no thermal destrnction), and ajv = 1. Up to the critical valne for A ~ 0.025, the 
agreement is very convincing. As we increase A into the mixed phase region of the diagram, 
the kinetic approach nnderestimates the amonnt of volnme in the +-phase. This is not 
snrprising, since for these valnes of A the density of snbcritical bnbbles is high enongh that 
terms not inclnded in the eqnation become important, as I mentioned before. However, the 
lack of agreement for higher valnes of A is irrelevant, if we are interested in having a measnre 
of the smoothness of the backgronnd; clearly, the rise in 7 is sharp enongh that homogeneons 
nncleation shonld not be trnsted for A > 0.024 or so, as the fraction of volnme occnpied by 
the +-phase is already aronnd 30% of the total volnme. Snbcritical bnbbles give a simple and 
qnantitatively accnrate pictnre of the degree of inhomogeneity of the backgronnd, offering a 
gnideline as to when homogeneons nncleation theory can be applied with conhdence. 


Figure 5: Comparison between results from numerical simulations (dots) and the result 
from the kinetic equation. 


B. Thermal Phase Mixing in the Standard Electroweak Model 


It is straightforward to compute 7 for the electroweak model. One can choose any 
effective potential, and simply plug the results into Eq. Here, I will work with the 
simplest version of the model, the unimproved 1-loop potential, written in Eq. Clearly, 
different effective potentials will give different values for 7 . So, if we used instead the 2-loop 
potential of Ref. p 6 |, we would obtain smaller values for 7 , as the two-loop potential gives 
rise to a stronger transition. The point here is not to get into the problems of computing a 
reliable effective potential, but to apply the method just developed to a simple situation. 

There have been several recent papers on applying snbcritical bubbles to different effec¬ 
tive potentials or exploring some of their properties |^. In my view, there are two important 
points which seem to be overlooked in most (but not all) of these works. First, the impor¬ 
tance of snbcritical bubbles is sensitive to the parameters of the model. If one hnds that 
snbcritical bubbles are not relevant for a given range of parameters, this does not mean that 
they are ruled out as unphysical. It simply means that they will be relevant for a different 
range (usually higher Higgs masses if one is interested in the electroweak phase transition); 
hot systems fluctuate, even if our description of the physics through a perturbative effective 
potential may be inaccurate. For the particular case of the electroweak transition, we ea- 
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gerly await for lattice gauge simulations which will probe the regime of higher Higgs masses. 
Second, approaches that use a saddle-point evaluation of the partition function to compute, 
say, the dispersion of the held around the background, will miss the important contributions 
coming from large amplitude huctuations. Thus, they will underestimate the amount of 
mixing. This is an advantage of the kinetic approach presented here. 

Using the parameters for the electroweak model, we can express 7 in terms of the Higgs 
mass, as shown in Fig. 6 . I took i?min = ^(^c), and 0c = 0+(Tc). Again, the results are 
sensitive to the choice of minimum radius and amplitude, but these values are certainly 
conservative. It is clear from the Figure that taking u = 1, 7 changes sharply from 0 to 
0.5 for 60 < Mniggs < 70, the window where lies the present lower bound on the Higgs 
mass. Thus, these results appear to rule out nucleation for the unimproved 1-loop effective 
potentiaQ. Mechanisms for baryogenesis based on homogeneous nucleation must take the 
mixing of the background due to subcritical fluctuations into account. 


Figure 6: Thermal phase mixing for the minimal electroweak model as a function of 
the Higgs mass. 


IV. CONCLUDING REMARKS 


In this talk I presented the results of both numerical simulations and analytical modelling 
of phase mixing induced by thermal fluctuations. Although I stressed mostly the importance 
of these results in the context of homogeneous nucleation in hnite temperature held theory, 
the ideas presented here can be adapted to a variety of different situations of interest also 
in laboratory applications, as long as the system in question can be described by a coarse¬ 
grained Ginzburg-Landau free energy with a non-conserved order parameter. For example, 
together with Andrew Heckler, a method to compute the inhuence of subcritical bubbles 
on the nucleation rate is presently being developed, which, we believe, is quite general . 


**There is a volume factor which suppresses 7 not included here; simply, not all of the subcritical 
bubble is in the -|—phase, just a fraction of it. However, this suppression is compensated by taking 
the more realistic choice 0c = 0 max 5 as opposed to the conservative choice 0c = 0+ we used. 
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The basic idea is to include the available free energy in the gas of subcritical bubbles in 
the computation of the effective nucleation barrier. In practice, this is done by obtaining a 
new coarse-grained effective action which sums over fluctuations of subcritical size, in the 
spirit of the renormalization group. Comparing the results from our method to numerical 
simulations of nucleation in 2 dimensions (Ref. [|^) shows how subcritical fluctuations can 
account for a dramatic drop in the effective nucleation barrier which is observed close to 
degeneracy. 

I hope it was clear from this talk that subcritical bubbles can provide a clear and quanti¬ 
tatively accurate description of the degree of homogeneity of a given background, metastable 
or not. It would be interesting to generalize the ideas presented here to quantum fluctuations 
as well, since we should expect them to be dominant at very low temperatures. Naively, 
we could think of quantum subcritical bubbles as nonperturbative violations of energy con¬ 
servation, with a probability distribution suppressed by their Euclidean action. If present, 
these fluctuations could have many interesting consequences, from inflationary scenarios to 
nucleation of topological defects in systems with nontrivial vacuum topology. 


Acknowledgements 


I am grateful to my collaborators Julian Borrill, Andrew Heckler and Rocky Kolb, as well 
as to M. Alford, E. Copeland, H.-R. Miiller, and R. Ramos for the many long discussions on 
bubbles and phase transitions. I am also grateful to the Nasa/Fermilab Astrophysics Center 
for their kind hospitality where some of the work presented here was developed. 

At Fermilab I acknowledge partial support from NASA Grant NAG 5-2788 and DOE. 
This work was partially supported at Dartmouth by the National Science Foundation 
through a Presidential Faculty Fellows Award no. PHY-9453431 and by Grant no. PHYS- 
9204726, and by a NASA Grant NAGW-4270. 

Finally, I thank Hector de Vega and Norma Sanchez for their warm hospitality during 
this workshop. 


References 

1. For a review see E. W. Kolb and M. S. Turner,. The Early Universe^ (Addison-Wesley, 
Redwood, GA, 1990). 

2. A. G. Gohen, D. B. Kaplan, and A. E. Nelson, Annu. Rev. Nucl. Part. Sci. 43, 27 
(1993); A. Dolgov, Rhys. Rep. 222, 311 (1992). 

3. See, however, G. R. Farrar and M. E. Shaposhnikov, Rhys. Rev. Lett. 70, 2833 (1993); 
71, 210(E) (1993); Rhys. Rev. D50, 774 (1994). 


18 


4. B. Liu, L. McLerran, and N. Turok, Phys. Rev. D46, 2668 (1992); A. G. Cohen, D. B. 
Kaplan, and A. E. Nelson, Phys. Lett. B336, 41 (1994); M. B. Gavela, P. Hernandez, 
J. Orloff, and O. Pene, Mod. Phys. Lett. A9, 795 (1994); P. Huet and E. Sather, Phys. 
Rev. D51, 379 (1994). 

5. R. Becker and W. Boring, Ann. Phys. 24, 719 (1935). 

6 . J. W. Cahn and J. E. Hilliard, J. Chem. Phys. 31, 688 (1959). 

7. J. S. Langer, Ann. Phys. (NY) 41, 108 (1967); ibid. 54, 258 (1969). 


8 . J. D. Gunton, M. San Miguel and P. S. Sahni, in Phase Transitions and Critical 
Phenomena, Vol. 8, Ed. C. Bomb and J. L. Lebowitz (Academic Press, London, 
1983). 

9. M. B. Voloshin, 1. Yu. Kobzarev, and L. B. Okun’, Yad. Fiz. 20, 1229 (1974) [Sov. J. 
Nucl. Phys. 20, 644 (1975) ]. 

10. S. Coleman, Phys. Rev. D15, 2929 (1977); C. Callan and S. Coleman, Phys. Rev. D16, 
1762 (1977). 

11. A. B. Linde, Phys. Lett. 70B, 306 (1977); Nucl. Phys. B216, 421 (1983); [Erratum: 
B223, 544 (1983)]. 

12. M. Gleiser, G. Marques, and R. Ramos, Phys. Rev. D48, 1571 (1993); B. E. Brahm 
and C. Lee, Phys. Rev. D49, 4094 (1994); B. Boyanovsky, B. E. Brahm, R. Holman, 
and B.-S. Lee, Nucl. Phys. B441, 609 (1995). 

13. P. Arnold and L. McLerran, Phys. Rev. D36, 581 (1987); ibid. D37, 1020 (1988). 

14. M. Gleiser, Phys. Rev. Lett. 73, 3495 (1994). 

15. J. Borrill and M. Gleiser, Phys. Rev. D51, 4111 (1995). 


16. M. Gleiser and R. Ramos, Phys. Rev. D50, 2441 (1994); B. L. Hu, J. P. Paz and Y. 
Zhang, in The Origin of Structure in the Universe, Ed. E. Gunzig and P. Nardone 
(Kluwer Acad. Publ. 1993); B. Lee and B. Boyanovsky, Nucl. Phys. B406, 631 (1993); 
S. Habib, in Stochastic Processes in Astrophysics, Proc. Eighth Annual Workshop in 
Nonlinear Astronomy (1993). 


17. See work by M. Gleiser and R. Ramos, in Ref. [|T^ 


18. M. Alford and M. Gleiser, Phys. Rev. D48, 2838 (1993). 


19. CBF Collaboration, F. Abe et al, Phys. Rev. Lett. 73, 225 (1994); Phys. Rev. D50, 
2966 (1994). 


20. See for example, N. Goldenfeld, Lectures on Phase Transitions and the Renormalization 
Group, Frontiers in Physics, Vol. 85, (Addison-Wesley, 1992). 


19 



21. M. Gleiser, E. W. Kolb, and R. Watkins, Nucl. Phys. B364, 411 (1991); G. Gelmini 
and M. Gleiser, Nucl. Phys. B419, 129 (1994); M. Gleiser and E. W. Kolb, Phys. Rev. 
Lett. 69, 1304 (1992); N. Tetradis, Z. Phys. C57, 331 (1993). 


22. G. Gelmini and M. Gleiser in Ref. 


0 . 


23. E.D. Siebert and G.M. Knobler, Phys. Rev. Lett., 52, 1133 (1984); J.S. Langer and 
A.J. Schwartz, Phys. Rev. A21, 948 (1980); A. Leggett in Helium Three, ed. by W.P. 
Halperin and L.P. Pitaevskii, (North-Holland, New York, 1990); for an (ontdated) 
review of the sitnation in the early eighties see Ref. |^. 


24. M. Gleiser, Phys. Rev. D49, 2978 (1994); E.J. Gopeland, M. Gleiser, and H.-R. Muller, 
DART-HEP-95/01, in press, Phys. Rev. D. 


25. M. Gleiser, A. Heckler, and E.W. Kolb, in progress. 


26. Z. Fodor and A. Hebecker, Nucl. Phys. B432, 127 (1994). 


27. M. Hindmarsh and R. Rivers, Nucl. Phys. B417, 506 (1994); T. Shiromizu, 
M. Morikawa, and J. Yokoyama, |hep-ph/9501312| ; K. Enqvist and 1. Vilja, [hep 


ph/9410224| ; J. Sirkka and 1. Vilja, |hep-ph/9503483| ; 1. Vilja, |hep-ph/9410224| ; K. 


Enqvist, A. Riotto, and 1. Vilja, Piep-ph/9505341|; F. Illuminati and A. Riotto, jhep- 


ph/9506419| ; L.M.A. Bettencourt, jhep-ph/ 9505440 


28. M. Gleiser and A. Heckler, in progress. 


20 


























TABLE 1 


A 

teq 

(7 

0.015 

25.0 

1.0 

0.020 

45.0 

1.0 

0.022 

60.0 

1.0 

0.024 

110.0 

0.80 

0.026 

220.0 

0.50 

0.028 

120.0 

0.75 

0.030 

100.0 

0.80 

0.035 

55.0 

0.80 

0.040 

40.0 

0.90 
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